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Abstract An expression for the square of the spin operator expectation value, is obtained for a 
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system. The optimization of the collinearity axis is also considered. 
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1 Introduction 

Particle-independent models based on single Slater determinant wave functions, have enjoyed consider¬ 
able interest in quantum chemistry, since the pioneering works of Hartree, Slater and Fock [1-3]. 

When a quantum system is described by a spin-free Hamiltonian, which obviously commutes with 
the spin operators Sz and a spin-symmetry respectfnl way of using the Hartree-Fock method consists 
in: 

1) Using spin-orbitals of pure a- or /3-spin, so that the HF optimized Slater determinant is an eigen¬ 
function of Sz', 

2) Imposing the spin-equivalence restriction [4], which means that paired a- and /3-spin-orbitals are 
formed from the same set of linearly independent orbitals. We have proved mathematically [5-7] that 
this additional constraint is a necessary and sufficient condition to insure that a Slater determinantal 
wave function is an eigenfunction of the spin operator S'^. In other words, we have shown that relaxing 
the 5^-symmetry constraint exactly amounts to allow different “paired orbitals”, in the sense of Refs. 
[8,9], to have different spins. This equivalence enabled us to characterize the variational space explored 
by the restricted open-shell Hartree-Fock (ROHF) method [10], which precisely consists in optimizing 
a Slater determinant subject to constraints 1) and 2) (plus spatial-symmetry constraints if any) [II]. 
The equivalence was also discovered independently [12] by optimizing a Slater determinant with a La¬ 
grange multiplier, enforcing (S^) to be arbitrarily close to the ROHF value, instead of applying the 
spin-eqnivalence restriction. Not surprisingly, the determinant was approaching the ROHF solution. 

Similar to the spin-free case is the “complex-free” one: When a quantum system is described by a 
real Hamiltonian, which obviously commutes with complex conjugation, one can restrict oneself to the 
calculation of real eigenfunctions. Then, it is also possible to employ only real spin-orbitals to construct 
the HF Slater determinant [13]. (However, difficulties may occur when the symmetry group of the 
molecule cannot be represented over real numbers, and nonetheless, one wishes the spin-orbitals to be 
adapted to spatial-symmetry). 

However, it has been proposed by various authors to relax some or all of the above-mentioned 
constraints, to gain variational freedom. For example, the different orbitals for different spins method 
(DODS) of Refs. [14,15], (which is usually just called “unrestricted Hartree-Fock” (UHF), but in this 
paper we use “DODS” to avoid confusions), relaxes the spin-equivalence restriction, hence the HF 
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solution is no longer an eigenfunction of S'^. Other authors [16-18] have advocated the use of general 
spin-orbitals, mixing a-spin and /3-spin parts, in conjunction with the use of projectors [19]. 

Along the same line of thought, the use of complex spin-orbitals has been proposed [17,20] to increase 
variational freedom in the case of real Hamiltonian. Prat and Lefebvre went a step further with so-called 
“hypercomplex” spin-orbitals to construct Slater determinants of arbitrary accuracy [21]. However, the 
coefficients of their spin-orbitals were elements of a Clifford algebra of dimension 2^”, that was not a 
normed division algebra, also known as Cayley algebra, for arbitrary values of n. This was unfortunate, 
since such a structure appears to be a minimal requirement for a quantum formalism, if, for example, 
Born’s interpretation of the wave function is to hold firmly. For n = 1, the Clifford algebra of Prat et al. 
was actually the non-commutative field of Quaternions, therefore, a fortiori, a normed division algebra. 
The only larger normed division algebra is the Octonion algebra. It is a Clifford algebra of dimension 
8, which has also been proposed in a quantum mechanical context [22], but this algebra is neither 
commutative nor associative. The lack of these properties rises difficulties for its use for multipartite 
quantum systems, nevertheless these difficulties can be overcome by keeping the product of octonion 
coefficients in the form of a tensor product. So, octonion-unrestricted HF appears to be the largest 
Clifford algebra-unrestricted single determinantal method that can be considered in the spirit of Prat 
and Lefebvre’s proposal. However, octonions seem incompatible with the desirable requirement that 
the algebra of quantum observables be what is now called a formally real Jordan algebra [23] acting 
on a vector space of arbitrarily large dimension. Octonions are also ruled out by the requirement of 
orthomodularity in infinite dimension according to Soler’s theorem [24,25], which restricts quantum 
Hilbert spaces to be real, complex or at most quaternionic. 

The first HF molecular calculations with general complex spin-orbitals, without projecting out the 
symmetry-breaking part of the wave function, are maybe those of Ref. [26]. It was found on the BH 
molecule around its equilibrium geometry that the general complex Hartree-Fock (GCHF) energy was 
indeed lower than the DODS one, which itself was lower than the restricted Hartree-Fock (RHF) solu¬ 
tion. So necessarily, the corresponding GCHF wave functions had S^-spin contamination and S'z-spin 
contamination, that is to say, the expectation values of these operators were different from 0, the value 
expected for a singlet ground state. (It is not clear whether complex numbers were used for this molecule, 
but the authors did mention that they performed complex calculations for 2-electron systems.) 

Relaxing the “Sz-constraint” hence the “collinearity constraint”, becomes perfectly legitimate when 
hyperfine or spin-orbit couplings are considered, since the operator Sz no longer commutes with the 
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Hamiltonian. As a matter of fact, real physical systems do exhibit either light [27,28] or strong [29, 
30] noncollinearity of their spin densities. Similarly, the use of complex spin-orbitals is natnral, when 
considering relativistic corrections resulting in a complex Hamiltonian operator. So, in such a context, 
one should use no less than general complex spin-orbitals in HF calculations [31]. The “spin-same-orbit” 
coupling term used in these calculations does not commute with the S'^-operator. Therefore, one cannot 
strictly speak of “S'^-spin contamination” in relativistic GCHF wave functions. However, calculating the 
expectation value of , a bona fide quantum observable, can still provide valuable physical information 
about the system. 

A general expression for the expectation value of has been obtained in the DODS case [8], and 
has served as a measure of S'^-spin contamination. However, as far as we are aware, no such formula has 
been published in the case of a GCHF wave function. This gap will be filled in the next section. 

Studying departure from collinearity is more difficult because of arbitrariness in the quantification 
axis. One possible way to overcome the difficulty would be to apply an external magnetic field to fix 
the z-axis but small enough not to perturb the GCHF solution. However, an elegant alternative has 
been proposed recently by Small et al. [32]. It is based on studying the lowest eigenvalue of a (3 x 3)- 
matrix built from expectation values of spin operator components and their products. In the GCHF 
case, the authors provided the expressions required to compute the matrix elements in a compact 
form. In the third section, we give a more extended formula in terms of molecular orbital overlap 
matrix elements. We also illustrate the connections between spin contamination, noncollinearity and its 
correlative: “perpendicularity” on the H 2 0^ cation example. We sum up our conclusions in the last 
section. 


2 Spin contamination in GCHF 

A General Complex Hartree Fock (GCHF) wave function 


^GCHF = (j)l A ■■■ A (j)N^ (1) 

is the antisymmetrized product (or wedge product, denoted by A) of orthonormal spinorbitals, or “two- 
component spinors”. 
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(3) 

where the scalar product (-I-) means integration over space variables and summation (i.e. taking the 
trace) over spin variables: {4>i\(j)j) = {4>ia\4>ja) + {4‘ip\4>j0)^ (where the same bracket symbol is used for 
the scalar product between orbital parts). We define the “number of a-spin electrons” (respectively 

Ne 

“number of /3-spin electrons”) as Na := X) ('/’fa(respectively, := ^ {4>i0\4'i0))■ It is the ex- 

i —1 2—1 

pectation value of the projection operator on the a- (respectively /3-) one-electron Hilbert subspace 
(more rigorously speaking, the operator induced onto the n-electron Hilbert space by this one-electron 
projection operator). Note that these two numbers need not be integer numbers, however their sum is 
an integer: Na + Np = Ng. 

Let us work out the expectation value of the spin operator. 


on a general GCHF wave function. 
The action of Sz is given by, 


S^ = S^, + -{S+S- + S-S+), 


1 ' 

Sz'l’GCHF = a Ac 


where, 


^hcHF — </>i A • • • A 4>i-i A A ipi+i A • • • A 


Note that, 


So, the expectation value of Sz is 




-| ^ -1 ^ .j ^ 

{^GCHf\Sz\^GCHf} = a (^GCgF « X! 1^*“^ “ {<t>ip\<t)i0)) = - -■ 


( 9 ) 
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and that of Si: 

{'^gchf\S1\(1>gchf) = {Sz*^gchf\Sz‘^gghf) 


Ne 


— Z (^gchfI^gchf) 


i,j=i 


Ne 

E 

i=l 


Ne 


= I I E{’^gghf\’^gghf) + E {^gchfI^gchf) 

i,j = l 


Ne 


Ne 


= iE (</'*!</'*) +E(-i)K</'*l</'i)P + 


i=l 


3 = 1 

i/i 


Ne 


— 4 Ne S- E ( ^)\{4’ia\4'ja) 'S {{4'ia\4’ia) {4'il3\4’if})) {{4’ja\4'ja) {4’jl3\4'jf})) 


i J = 1 

Ne 


— 4 1 -^e + E {{4’ia\4'ia) {4'i/3\4’if})) {{4'ja\4’ja) (^^i/3 )) \{4’ia\4'ja) (0i/3 I'/’j/S) I 

V *’^=1 y 


( 10 ) 


This equation reduces to in the case of a DODS wave function. So, the second term on 

the right-hand side (rhs), which is {{<1>gchf\S1\<1’gchf) — {‘1>gchf\Sz\‘1>gchf)^) is directly related to 
relaxation of the S'z-constraint and will be called the “^-noncollinearity” contribution. Note, however, 
that for a GCHF wave function, the first term on the rhs does not necessarily correspond to an eigenvalue 
of Si, according to the definition of and Np. 

The action of 5'+ is given by, 

S+<Pc 


Ne 


'^GCHF — ^GCHFy 
i=l 


where. 


and, 


^GCffF — (/>i A • • • A A (f>i A (j)i+i A ■ ■ ■ A (j)Ne 

+ 4>iP 

0 


( 11 ) 


( 12 ) 


(13) 


Similarly, the action of S is given by. 


Ne 


S ^GCHF = y^.^GCHF^ 

i=l 


( 14 ) 
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where, 


^hcHF — </>i A • • • A A A (j)i+i A • • • A 


(15) 


and, 


0 


(16) 


So, the expectation value of S S~^ is. 


{'1^gchf\S~ S~^\<Pgchf) = {S~^<Pgchf\S~^<Pgchf) 

Ne ,, 

= Z) {^gchfI^gchf) 

i,i=l 

iVe , . , . 

= Z (^GCffFl^GCffF) + Z (^GCFfI^GCFf) 

i=l 

iAi 

Ne ( , , Ne 

= Z + Z(-i)l(^#i)P + (^il<)'*)('?^'il^j) 

i=l \ 3=1 

\ i/i 

Ne I Ne 

~ Z I {4’i/3\'Pif3) A- '^{~^)\{<Pi/3\4’ja)\ + | ^ice) | ^j/3 ) 

i=l \ 3=1 

\ 

Ne 

— -A^/3 + Z (^^i/3 l^ia) l^^i/s) ■ 

i,i=l 


(17) 


Similarly, the expectation value of S'^S is. 


{'1^gghf\S~^S |<?gcff) = {S <Pgchf\S <I>gchf) 

Ne I ,, Ne 

= Z {(|>^\(|>^) + Z(-i)l(^il«^'3)P + 

i=l \ 3=1 

\ JAi 


Ne 


Ne 


= z (0za|0za)H“ X] ( 1) I I {4^ia.\4^i0) {4^j^\4^ja) 


3 = 1 
3^^ 


Ne 


— “1“ {4^ict\4^il3') {4^j 0\4^ja'} (^zct | ) (^j/3 | ^zct) ■ 


( 18 ) 
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Using Eq.(4) and putting together Eqs.(lO), (17) and (18), one obtains the expectation value of 3“^, 


Ne 

+ X/ {4^ia\<pip){<j>j0\4^ja) ~ {<l>ia\4’j/3){4’j/3\4’ia)- (19) 

jj = l 

The expression reduces to the known formula in the case of a DODS wave function. Assuming, without 
loss of generality, that > N^, we rewrite Eq. (19) as, 

{^GCHf\S‘^\^GCHf) = - (^-Y -^ 4 ^ I' 

( \ Ne 

+ lV/3 - X) + I E 1 ) P ■ (20) 

\ *J = 1 / i=l 


In this formula we identify four contributions: The first term is formally identical to the ROHE expression 
also found in the DODS case. However, care must be taken that it is actually different, because the 
numbers of a- and /3-electrons are not good quantum numbers in the GCHF case. The second term 
on the first line is the “z-noncollinearity” contribution. The third term on the second line is formally 
analogous to the “spin contamination” of a DODS wave function as defined in [7,8]. Finally, the last 
term on the second line is the square of the expectation value of the lowering or raising operator: 

Ne 

= \{^gchf\S^\^gchf)\^ = \{^gghf\S~\^gghf)\^ ■ (21) 

A non zero contribution of this term can only arise from the release of the Sz-constraint, which allows 
for the a- and /3-components of a given, general spin-orbital to be both non zero. But it originates 
from S~^S~ and S~S~^, and is maximal when (pip = exp(z0)(/)ia for all i, that is to say when the pi's 
are eigenfunctions of cosOSx + sinOSy for some angle 9. It is related to the emergence of a non-zero 
spin density in the x, y-plane, correlatively to the loss of z-collinearity. We tentatively call this term the 
“x, //-perpendicularity” contribution. 

The present formulas have been implemented in the code TONTO [33] and applied in a recent 
article (third column of Tab. 4 in [34]). Let us discuss further the different contributions to for a 
7720+ GCHF calculation similar to that reported in [34]. The z-quantification axis was the axis per¬ 
pendicular to the plane of the molecule. The results, see Tab. 1, shows that the main contribution to 
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{'^GCHf\S'^\'^gchf) beside the reference expression (first term in eq.(20)) is the so-called spin contami¬ 
nation contribution (we set h = 1 throughout the paper). The x, j/-perpendicularity and z-noncollinearity 
contributions are of the same order of magnitude and more than one order of magnitude smaller. Added 
~ ■^) ^ + l)’ almost make up the reference value of -1-0.75. So, the spin con¬ 

tamination value of 0.007033 amounts almost exactly to the difference between the exact expectation 
value {'1>gghf\S^Wgghf) and this reference “ROHF value”. This demonstrates that, in Table 4 of 
Ref. [34], the equality of the entries in column 2 (reference ROHF value plus spin contamination term) 
and column 4 (our {'1’gghf\S‘^\'1’gghf) value) does not imply no noncollinearity. In contrast, if for a 
given line of the table, these two quantities differ, then necessarily there will be some noncollinearity 
in the corresponding GCHF wave function. This can be shown reductio ad absurdum. Suppose that the 
z-collinearity constraint is fullfiled, then and Np will be good quantum numbers and the first term 
in our expression of {d>GCHF\S^\^GCHF) will be equal to the ROHF reference value. The spin contami¬ 
nation contribution being included in both quantities, the difference between them must arise from the 
X, y-perpendicularity and 2 -noncollinearity contributions. One at least of the contributions arising from 
the release of the collinearity constraint must be non zero, hence a contradiction. This hints that the 
following systems Cl, HCC, Fe, Cu, Cu^~^ and [OsCl 5 {Hpz)]~ reported in Table 4 of Ref.[34] would 
present stronger non collinearity than H 20 '^. 


3 Collinearity in GCHF 

In the previous section, we have encountered a 2 ;-(non)collinearity measure, colz '■= {{8^) — (*S'z)^). This 

^ Uj,\ 

quantity can be generalized to an arbitrary quantization direction defined by a unit vector u = u, 

[uzj 

of the unit sphere 5^ of by replacing S'z by u • S = ^ Uf^S^. Then, u-(non)collinearity is 


measured by: 


co/(u) ^ - {Sfj,){S„)). 


( 22 ) 


Small et al. [32] defined a (non)collinearity measure by: 


col := min col(u), 


(23) 
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which corresponds to the lowest eigenvalue of the matrix A whose elements are given by, 




(24) 


where is the real part of z. The associated eigenvector gives the optimal collinearity direction. 

Setting + (/>), y(j) = ^{4> — (j)) and ^(j) = cj), we have in this notation, 

N 

€ {x,y,z} A^^ (25) 

i,i=l 


where is the Kronecker symbol. 


Returning to the H20'^ example and applying these formulae, we obtain 


A = 


^+0.253128 +0.000145 -0.009774^ 
+0.000145 +0.253451 +0.003745 
j^-0.009774 +0.003745 +0.000461y 


The diagonalization of the ^-matrix gives the optimal collinear direction: 


uo‘ = (+0.0385908,-0.014789,+0.999146), 


(26) 


(27) 


which is only slightly tilted with respect to the z-direction, and the system is quasi-collinear in this 
direction since col = 0.000028 is very close to zero. This shows that the noncollinearity contribution 
to {'1'gchf\S^\<Pgchf) could be further reduced by more than one order of magnitude by selecting 
the optimal quantization axis corresponding to uo instead of the spatial z-axis. The perpendicularity 
contribution would decrease accordingly. 


4 Conclusion 

We have decomposed the expectation value of the spin operator into (i) a term formally identical to 
its expression for a ROHF reference wave function, (ii) a term called “spin contamination” because it is 
formally analogous to that derived by Amos and Hall [8] for DODS wave functions, (iii) a noncollinear 
contribution which can be minimized by following a procedure recently introduced [32], (iv) a term 
called the “perpendicularity contribution” which arises from the release of the z-collinearity constraint 
but which should rather be regarded as arising from the release of the “nonperpendicularity constraint” 
on the spin-density. The collinearity and nonperpendicularity constraints are correlatives. 
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We have evaluated these four different contributions for a GCHF calculation on the H20~^ cation. 
Note that we used the lOTC relativistic Hamiltonian [36] so that the term “spin contamination” is 
not really appropriate in this context, departure from the ROHF reference value being legitimate. 
However, the so-called spin contamination contribution has been found to dominate the noncollinearity 
and perpendicularity ones. This could be made even more so, by tilting the quantification axis to the 
optimal collinearity direction. 
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+4.999546 

Np 

+4.000454 


+0.749091 

2 -noncollinearity 

+0.000461 

X, y-nonperpendicularity 

+0.000427 

spin contamination 

+0.007033 

(S2> 

+0.757013 


Table 1 Expectation value of (S'^) and related quantities for an H20"*" GCHF optimized wave function. The geometry 
parameters were roH = 0.99192A, HOH = 101.411de^. The basis set consisted of the primitives Gaussian functions 
left uncontracted of Dunning’s cc-pVDZ hydrogen and oxygen basis sets [35]. The infinite-order two-component (lOTC) 
relativistic Hamiltonian of Barysz and Sadlej [36] was employed. 





